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ABSTRACT 

Recent observations show that the temporal decay of the R-band afterglow from GRB 990123 steepened 
about 2.5 days after the burst. We here propose a possible explanation for such a steepening: a shock 
expanding in a dense medium has undergone the transition from a relativistic phase to a nonrelativistic 
phase. We find that this model is consistent with the observations if the medium density is about 
3 X 10^ cm""^. By fitting our model to the observed optical and X-ray afterglow quantitatively, we 
further infer the electron and magnetic energy fractions of the shocked medium and find that the two 
parameters are about 0.1 and 2 x 10~^ respectively. The former parameter is near the equipartition value 
while the latter is about six orders of magnitude smaller than inferred from the GRB 970508 afterglow. 
We also discuss possibilities that the dense medium can be produced. 

Subject headings: gamma rays: bursts - shock waves 



1. INTRODUCTION 

The gamma-ray burst (GRB) 990123 was an extraordi- 
' nary event. It was the brightest burst yet detected with 
, the Wide Field Camera on the BeppoSAX satellite (Fe- 
' roci et al. 1999), and had a total gamma-ray fluence of 
\ ~ 5 X 10""* ergcm"^, which is in the top 0.3% of all bursts. 
■ It was the first burst to be simultaneous detected in the 
optical band. Optical emission with peak magnitude of 
^ 9 was discovered by the Robotic Optical Transient 
Search Experiment (ROTSE) during the burst and was 
' found to have rapidly faded down immediately after the 
. gamma-ray emission (Akerlof et al. 1999). The detection 
' of the redshift showed that the burst appears at z > 1.6 
; (Andersen et al. 1999; Kulkarni et al. 1999a). This im- 
. plies that if the GRB emission was directed isotropically, 
the inferred energy release is > 1.6 x 10^^ ergs (Kulkarni 
et al. 1999a; Briggs et al. 1999). 

The burst's afterglow was detected and monitored at 
' X-ray, optical and radio bands. It was the brightest 
of all GRB X-ray afterglows observed until now. The 
BeppoSAX detected the flux of the afterglow at 2-10 
keV six hours after the gamma-ray trigger to be 1.1 x 
10~^^ ergcm^^ s^^ and the subsequent temporal decay in- 
dex to be ax = -1A4 ± 0.07 (Heise et al. 1999a, b). The 
R-band optical afterglow about 3.5 hours after the burst 
showed a power-law decay with index am = — 1.1 ± 0.03 
(Kulkarni et al. 1999a; Castro-Tirado et al. 1999; Fruchter 
et al. 1999). This law continued until about 2.04 ± 0.46 
days after the burst. Then the optical emission be- 
gan to decline based on another power law with index 
a2R. = -1.65±0.06 (Kulkarni et al. 1999a) or -1.75±0.11 
(Castro-Tirado et al. 1999) or -1.8 (Fruchter et al. 1999). 
In addition, a radio flare was also detected about 1 day af- 
ter the burst (Kulkarni et al. 1999b; Galama et al. 1999). 

A scenario has been proposed to explain these observa- 
tions. If the burst is assumed to be produced from a jet, 



the steepening of the late optical afterglow decay is due to 
the possibility that this jet has undergone the transition 
from a spherical-like phase to a sideways-expansion phase 
(Rhoads 1997, 1999; Kulkarni et al. 1999a; Fruchter et 
al. 1999; Sari, Piran & Halpern 1999) or that we have 
observed the edge of the jet (Panaitescu & Meszaros 1998; 
Meszaros & Rees 1999). 

In this Letter we propose another possible scenario, in 
which the steepening of the late optical afterglow decay 
is due to the shock which has evolved from a relativis- 
tic phase to a nonrelativistic phase in a dense medium. 
According to the standard afterglow shock model (for a 
review see Piran 1998), the afterglow is produced by syn- 
chrotron radiation or inverse Compton scattering in the 
external forward wave (blast wave) of the GRB fireball 
expanding in a homogeneous medium. The external re- 
verse shock of the fireball may lead to a prompt optical 
flash (Sari & Piran 1999). As more and more ambient 
matter is swept up, the forward shock gradually decel- 
erates and eventually enters a nonrelativistic phase. In 
the meantime, the emission from such a shock fades down, 
dominating at the beginning in X-rays and progressively at 
optical to radio energy band. There are two limiting cases 
(adiabatic and highly radiative) for the hydrodynamical 
evolution of the shock. These cases have been well stud- 
ied both analytically (e.g., Meszaros & Rees 1997; Wijers, 
Rees k Meszaros 1997; Waxman 1997a, b; Reichart 1997; 
Sari 1997; Vietri 1997; Katz & Piran 1997; Meszaros, Rees 
& Wijers 1998; Dai & Lu 1998a; Sari, Piran & Narayan 
1998; etc) and numerically (e.g., Panaitescu, Meszaros & 
Rees 1998; Huang et al. 1998; Huang, Dai & Lu 1998). 
A partially radiative (intermediate) case has been investi- 
gated (Chiang & Dermer 1998; Cohen, Piran & Sari 1998; 
Dai, Huang & Lu 1999). Here we only consider the limit- 
ing cases. In the highly radiative model, since all shock- 
heated electrons cool faster than the age of the shock, the 
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optical afterglow should have the same temporal decay 
index as the X-ray afterglow (Sari ct al. 1998), incom- 
patible with the observations (Kulkarni et al. 1999a). In 
the adiabatic model, however, the difference in the decay 
index between optical and X-ray afterglows is found to be 
likely 1/4, which is consistent with the observational re- 
sult Aa = am — ax ~ 0.3. This implies that the shock 
producing the afterglow of GRB 990123 has evolved adia- 
batically. This is the starting point of our analysis. For an 
adiabatic shock, the time at which it enters a nonrclativis- 
tic phase oc n^^/"*, where n is the baryon number density 
of the medium. Therefore, this time for a shock expanding 
in a dense medium with density of n ~ 10^ cm~^ is two 
orders of magnitude smaller than that for a shock with the 
same energy in a thin medium with density of n ~ 1 cm^''. 
Furthermore, as given in Section 2, the afterglow at the 
nonrelativistic phase decays faster than at the relativistic 
phase. It is natural to expect that this effect can provide 
an explanation for the steepening feature of the afterglow 
from GRB 990123. 

Dense media have been discussed in the context of 
GRBs. First, Katz (1994) suggested collisions of relativis- 
tic nucleons with a dense cloud as an explanation of the 
delayed hard photons from GRB 940217. Second, to ex- 
plain the radio flare of GRB 990123, Shi & Gyuk (1999) 
speculated that a relativistic shock may have ploughed 
into a dense medium off the line of sight. Third, Piro 
et al. (1999) and Yoshida et al. (1999) have reported 
an iron emission line in the X-ray afterglow spectrum of 
GRB 970508 and GRB 970828 respectively. The observed 
line intensity requires a dense medium with a large iron 
mass concentrated in the vicinity of the burst (Lazzati, 
Campana & Ghisellini 1999). Finally, dense media (e.g., 
clouds or ejecta) may appear in the context of some en- 
ergy source models, e.g., failed supernovae (Woosley 1993), 
hypernovae (Paczynski 1998), supranovae (Vietri & Stella 
1998), phase transition of neutron stars to strange stars 
(Dai & Lu 1998b), baryon decay of neutron stars (Pen & 
Loeb 1998), etc. 

2. THE EVOLUTION OF A SHOCK IN A DENSE MEDIUM 

2.1. Relaiivistic Phase 

Now we consider an adiabatic relativistic shock expand- 
ing in a dense medium. The Blandford-McKee (1976) self 
similar solution gives the Lorentz factor of the shock. 
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where E = E54 x lO^^ergs is the total isotropic energy, 

1^5 = n/10^ cm~^, = iday x 1 day is the observer's time 
since the gamma-ray trigger, z is the the redshift of the 
source generating this shock, and mp is the proton mass. 

In analyzing the spectrum and light curve of synchrotron 
radiation from the shock, one needs to know two cru- 
cial frequencies: the synchrotron radiation peak frequency 
i^m) and the cooling frequency (t^c)- In the standard after- 
glow shock picture, the electrons heated by the shock are 
assumed to have a power-law distribution: dNe/dje oc 
for 7e > 7em) where 7e is the electron Lorentz factor and 



the minimum Lorentz factor 7em — GlOceT- The power- 
law index p « 2.56 by fitting the spectrum and light curve 
of the observed afterglow of GRB 990123 (see below). We 
further assume that eg and eb are ratios of the electron 
and magnetic energy densities to the thermal energy den- 
sity of the shocked medium respectively. Based on these 
assumptions, the synchrotron radiation peak frequency in 
the observer's frame can be written as 
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where es-Q = es/lO"^ and B' = (327re_B7^nmpC^)-^/^ 
is the internal magnetic field strength of the shocked 
medium. According to Sari et al. (1998), the cooling fre- 
quency, the frequency of electrons with Lorentz factor of 
7c that cool on the dynamical time of the shock, is given 

by 
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where ctt is the Thompson scattering cross section. From 
equations (2) and (3), Sari et al. (1998) have further de- 
fined two critical times, when the breaking frequencies Vm 
and cross the observed frequency v = x 10^^ Hz: 
= 8.6 X lO-^ef/'e^/iei^s'f [(1 + ^)/2.6]V3^-2/3days, 

and tc = 380e^^_gi;5"/n^^[(l + z)l2.&f^v^^ days. There- 
fore we see that for £54 ~ 1.6, e^ ~ 0.1, es.-e ~ 0.02, 
and ns ~ 30 inferred in the next section, the optical 
afterglow in several days after the burst should result 
from those slowly-cooling electrons and the X-ray after- 
glow from those fastly-cooling electrons. 

The observed synchrotron radiation peak flux can be 
obtained by 



N,jP'{l + z) 
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where Ne is the total number of swept-up electrons, = 
rrieC^aTB' /{2,e) is the radiated power per electron per 
unit frequency in the frame comoving with the shocked 
medium, and Dl = -0^,28 x 10^^ cm is the distance to the 
source. In the light of equations (2)- (4), one can easily 
find the spectrum and light curve of the afterglow, 
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where the low-frequency radiation component has not been 
considered (Sari et al. 1998). In the GRB 990123 case, we 
require Vm < f < for the optical afterglow and u > Vc 
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for the X-ray afterglow. Thus, the R-band afterglow de- 
cay index = 3(1 — p)/4 and the X-ray decay index 
ax = (2 — 3p)/4, which are well consistent with the obser- 
vational results aiR = 1.1 ± 0.03 and ax = -1.44 ± 0.07 
if 2.56. 

2.2. Nonrelativistic Phase 

As it sweeps up sufficient ambient matter, the shock will 
eventually go into a nonrelativistic phase. During such a 
phase, the shock's velocity v oc f^^^^, its radius r oc t^^, 
the internal field strength B' oc f^^^^ and the typical elec- 



this result with equations (2)- (5), we can also derive 
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Thus, we obtain the 



synchrotron peak frequency I'm oc '^^m^' ^e"^; the cool- 

Q 2 1/5 

ing frequency Uc oc B ~ oc and the peak flux 

Fjy^ oc Nf.Pl^ oc r^B' oc t^^. According to these scaling 
laws, we further derive the spectrum and light curve at the 
nonrelativistic stage: 
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From equation (6), we can see the R-band decay index 
an = (21 — 15p)/10 for radiation from slowly-cooling elec- 
trons or an = (4 — 3p)/2 for radiation from rapidly-cooling 
electrons. If p « 2.56, then an w —1.74 or —1.84, in excel- 
lent agreement with the observations in the time interval of 
2.5 days to 20 days after the burst (Kulkarni et al. 1999a; 
Fruchter et al. 1999; Castro-Tirado et al. 1999). 

3. CONSTRAINTS ON PARAMETERS 

In the above section, we show that as an adiabatic shock 
expands in a dense medium from an ultrarelativistic phase 
to a nonrelativistic phase, the decay of the radiation from 
such a shock will steepen. This effect may fit the observed 
steepening better than the alternative interpretation — jet 
sideways expansion. In the latter interpretation, the tem- 
poral decay of a late afterglow is very likely to be cx 
(Rhoads 1997, 1999: Sari et al. 1999).' We further analyze 
our effect and infer some parameters of the model. 

According to the analysis on the R-band light curve of 
the GRB 990123 afterglow (Kulkarni et al. 1999a; Fruchter 
et al. 1999; Castro-Tirado et al. 1999), the observed 
break occurred at = 2.04 ± 0.46 days. This implies 
7 ~ 1 at tday ~ 2.5. From equation (1), therefore, we find 
ns ~ 16£'54, where the redshift z = 1.6 has been used. 
We now continue to consider two observational results. 
First, on January 23.577 UT, the Palomar 60-inch tele- 
scope detected the R-band magnitude i?, = 18.65 ± 0.04, 
corresponding to the flux Fji ~ 100 /xJy at fday ~ 0.17 
(Kulkarni et al. 1999a). Considering this result in equa- 
tion (5) together with equations (2) and (4), we can derive 
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where the right number has been obtained by taking 
p fa 2.56 and Dl,28 ~ 3.7. Second, on January 24.65 
UT, the BeppoSAX observed the X-ray (2-10 keV) flux 
Fx ^ 5 X 10~^/xJy (Heise et al. 1999a, b). Combining 
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Since E54 ~ 1.6 (Briggs et al. 1999; Kulkarni et al. 1999a), 
the medium density ~ 30 and the solution of equations 
(7) and (8) is Ce ~ 0.1 and es-e ~ 0.02. Our inferred 
value of Eg is near the equipartition value, in agreement 
with the result of Wijers & Galama (1998) and Granot, 
Piran & Sari (1998), while our es is about six orders of 
magnitude smaller than the value inferred from the after- 
glow of GRB 970508. Of course, the field density for GRB 
971214 has been estimated to be less than 10~^ times the 
equipartition value (Wijers & Galama 1998). As suggested 
by Galama et al. (1999), such differences in field strength 
may refiect differences in energy flow from the central en- 
gine. 

4. DISCUSSION AND CONCLUSION 

In the above section, we find the medium density 
n ^ 3 X lO^cm"^ for our model to fit the observed 
optical and X-ray afterglow of GRB 990123. Now we 
show that even in the the presence of such a dense 
medium, the optical and X-ray radiations from the for- 
ward shock were neither self absorbed in the shocked 
medium nor scattered in the unshocked medium. First, 
the self-absorption frequency of the shocked medium is 
(Wijers & Galama 1998; Granot et al. 1998) ~ 
10^ GHz(ee/0.1)-i(es,-6/0.01)i/%^f (ns/lO)^/^ This 
estimate should be the upper limit because of the pres- 
ence of a possible low-energy electron population (Wax- 
man 1997b). Clearly, Va is much less than the optical fre- 
quency, implying that the self absorption in the shocked 
medium didn't affect the optical and X-ray afterglow. In 
fact, this estimate is valid only for Va < Pm- When 
t'a > i^m, t^a niust have decayed. As a result, the flux 
at 8.46 GHz first increased as i^^^ and then declined as 

^0^ ''^ for Ua < 8.46 GHz during the nonrelativistic phase. 
This might provide an explanation for the observed ra- 
dio fiare. Second, a photon emitted from the shock may 
be scattered by the electrons in the unshocked medium. 
The scattering optical depth r ^ axnU (where R is the 
typical radius of the medium). If the medium was dis- 
tributed isotropically and homogeneously and its mass 
M ^ IOM0 (the typical mass of a supernova ejecta), then 
r ~ 0.05(M/10MQ)i/3(n5/10)2/3 < 1. This implies that 
the afterglow from the shock was hardly affected by the 
medium. 

For other well-studied afterglows, e.g., GRB 970228 and 
GRB 970508, their ambient densities must be very low for 
three reasons: (i) In these bursts there was no observed 
break in the optical light curve as long as the afterglow 
could be observed (Fruchter et al. 1998; Zharikov et al. 
1998). (ii) The fluctuation appearing in the radio after- 
glow light curve of GRB 970508 requires the shock had 
been relativistic for several weeks (Waxman, Kulkarni & 
Frail 1998) . (iii) The analysis of the afterglow spectrum of 
GRB 970508 leads to a low ambient density n < 10 cm"^ 
(Wijers & Galama 1998; Granot et al. 1998). However, 
the observed iron emission line in the X-ray afterglow spec- 
trum of GRB 970508 indeed requires a dense medium with 
density ~ 10^ cm"^ (Lazzati et al. 1999). The only way to 
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reconcile a monthly lasting power-law afterglow with iron 
line emission is through a particular geometry, in which 
the line of sight is devoid of the dense medium. In con- 
trast to this idea, we suggest that for GRB 990123 a dense 
medium of n ~ 3 x 10^ cm~^ appears at least at the line 
of sight or perhaps isotropically. 

How was the dense medium produced? One possibility 
was a cloud and another possibility was an ejecta from the 
GRB site. There have been several source models (men- 
tioned in Introduction) in the literature which may lead 
to massive ejecta. Here we want to discuss one of them 
in detail. Timmes, Woosley & Weaver (1996) showed that 
Type II supernovae may produce a kind of neutron star 
with ~ 1.73Mq. If these massive neutron stars have very 
short periods at birth, they may subsequently convert into 
strange stars due to rapid loss of angular momenta (Cheng 
& Dai 1998), and perhaps the strange stars are differen- 
tially rotating (Dai & Lu 1998b). Even though this model 
is somewhat similar to the supranova model of Vietri & 
Stella (1998), resultant compact objects are strange stars 
in our model and black holes in the supranova model. We 
further discuss implications of our model. First, the model 
leads to low-mass loading matter because of thin baryonic 
crusts of the strange stars. Second, such stars result in 
GRBs with spiky light curves, being consistent with the 
analytical result from the observed data of GRB 990123 
(Fenimore, Ramirez- Ruiz & Wu 1999). The third advan- 
tage of this model is to be able to explain well the prop- 
erty of the early afterglow of GRB 970508 by considering 
energy injection from the central pulsar (Dai & Lu 1998b, 
c). Finally, a dense medium, the supernova ejecta, appears 
naturally. 



Our scenario proposed in this Letter requires a dense 
medium with density 3 x 10^ cm^"^ to explain the 
steepening in the temporal decay of the R-band afterglow 
about 2.5 days after GRB 990123. We also suggest that 
this medium could be a supernova/supranova/hypernova 
ejecta. Thus, if the mass of the medium is assumed 
to be M ~ IOMq, its radius can be estimated to be 
i? - 3 X 10i'^cm(M/10MQ)i/3(„5/io)-i/3, According to 
equation (1), we can integrate dr — 27^cc?i® and thus find 
that the postburst 2.5-day time in the observer's frame 
corresponds to about 20 days in the unshocked medium's 
frame. This implies that the radius at which the shock 
entered a nonrelativistic phase is about 5 x 10"'^^ cm. This 
radius is much less than that of the medium. Therefore, 
the medium discussed here was so wide and dense that 
the ultrarelativistic shock must have become nonrelativis- 
tic about 2.5 days after the burst. 

In summary, a simple explanation for the "steepening" 
observed in the temporal decay of the late R-band after- 
glow of GRB 990123 is that a shock expanding in a dense 
medium with density of 3 x 10^ cm"'^ has evolved from a 
relativistic phase to a nonrelativistic phase. We find that 
this scenario not only explains well the optical afterglow 
but also accounts for the observed X-ray afterglow quan- 
tatitively. 
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